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1 INTRODUCTION 



1.1 X-ray Supernovae 



' Supernovae (SNe) are stellar explosions which are lumi- 
I nous in various frequencies of electromagnetic wave. Nowa- 
, days, the detection rate of SNe in Optical is about few 
' hundreds per year, while number of SNe detection in X- 
rays are still rare which is less than 50 detections in to- 
1 tal. There are two kinds of X-ray emission mechanisms 
in SNe which are circumstellar-shock interaction or shock 
breakout. In shock breakout case, shock breakouts enve- 
lope of the progenitor which emits strong X-rays in soft 
band for a short period of time, says minutes. X-ray tran- 
sient XRO 0801 9 or S N 2008D located in NGC 2770 (see 
ISoderberg et all (|2008D ) is the only case discovered. In cir- 
cumstellar interaction, the X-ray emission is relatively sta- 
ble which can be last for few years. Most of the radia- 
tions come from the explosion shock wave pair (i.e. "for- 
ward shock" and "reverse shock"). When shock from the 
SN center interacts with circumstellar matter (GSM) de- 
posited by stellar wind from the progenitors/progenitors' 
companion star. X-ray photons are produced d ue to thermal, 
inverse Compton and synchrotro n emissions (jlmmler et al.l 
l2006l : IChevalier fc Franssonll2006l). Accordi ng to the circum- 
stellar interaction model by IChevalier fc Franssoru (|l994| ). 
cooling shell and forward/reverse shocks pair are formed as 
the ejecta collides with the CSM. In most cases, the forward 
shocks heat up the CSM to emit hard X-ray (> 10 keV) 
while the reverse shocks interact with the ejecta to emit soft 
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X-ray (< 5 keV). In the early time, the X-ray opaque cool 
dense shell between the shocks blocks the soft X-ray from 
the reverse shocked region. After several days or weeks, the 
expanding cool shell becomes optically thin which allows 
the pass through of the soft X-ray. Since there is a rela- 
tively high density in the reverse shocked region, soft X-ray 
is strong and therefore dominates the entire X-ray emission. 
One of the oldest X-ray supernova detected was SN 1970G 
by a Chandra deep image o f MlOl in 2005, 35 years after 
the explosion (|Triimperll2008h . The luminosity of this SN was 
~ lO^'^ergs"^ in 0.3-2 keV which is consisten t with a mass 
loss rate of 2.6 (fu,/10 kms~^)(Af0/ yr s"'^) (|lmmler et al.l 
[20051). Recent study indicated that there were e ven older X- 
ray S Ne: SN 1941C, SN 1959D and SN 1968D (jSoria et all 



l2008h with X -ray luminosity around 10"^^ "^^erg s m energy 
range 0.3-8 keV which were detected by Ctiandra. These old 
supernovae are important because they fill up the black box 
between old and young phases of a supernova remnant. 



In this work. X-ray telescope (XRT) on Swift was used 
for searching X-ray SN candidate from different evolution 
phases. As the Swift archival data is plentiful (over 30 ob- 
servations in different field everyday), it is a suitable place 
to do such X-ray SNe survey. Total 24 X-ray SNe have been 
detected in this survey. 3 of them are newly-discovered, i.e. 
SN 1986L, SN 20031x and SN 2007od. Scientific statistical 
results of this survey are presented in §2.3 and individual 
studies of some candidates are discussed in 53. 
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2 RESULTS AND DATA ANALYSIS 
2.1 Searching Method 

There are over 5000 SNe have been opticaUy discovered since 
1885 which were recorded in detail in the Center for As- 
trophysics (CfA) archive 0. Since the sample size of survey 
target is so large, an automated program was developed to 
screen all Swift archival data in photon counting mode in 
the period November 2004 to February 2011 and detect if 
there is any X-ray SNe in the field. Here gives a brief idea 
about how the program was working. 

Firstly, we developed a local Swift database by down- 
loading all old archival data from the Swift public archive 
and then using a time-based job scheduler cron to update 
fresh calibrated X-ray observations to the database at every 
midnight. Since the file size of a full data set (i.e. auxiliary 
files, UVOT data and etc) is large which is usually larger 
than 10 Mb per observation in size, so only X-ray event files 
which contain information of arrival photons were down- 
loaded to minimize time and computer resources required. 
All observations in the database were classified into groups 
according to their aim points. 

Afterward, A list which contains all detected SNe po- 
sitions was extracted and updated from CfA. The program 
then did positional matching between the the list and the 
database. Since the real field of view of Swift is 23'. 6 x 23'. 6, 
we set the critical separation to be half length of the diagonal 
r = 16'. 7. Once the angular separation between a observa- 
tion and a SN is less than the critical value, the existence of 
the corresponding SN in the field of view is confirmed. 

For observations that contain SN/SNe in their fields of 
view, full set of data were downloaded. Observations that 
belong to the same SN were merged by xselect to optimize 
the signal to noise ratio and a corresponding X-ray image 
in energy range 0.3 to 10 keV was extracted. All extracted 
images were then processed by wavdetect of ClAO to detect 
all X-ray sources in the field. We set the wavelet scales to be 
a series from 2 to 8 increasing by a factor of %/2 , (i.e. 2.0, 
2.828, 4.0, 5.657, 8.0). In order to detect those X-ray sources 
with low counts, the values of the signal and background 
thresholds were set to 5 x 10~^. In addition, exposure maps 
generated by xrtexpomap were used to correct the variation 
of effective area throughout the image. 

According to the pointing accuracy and the angular res- 
olution of Swift XRT, a box search of size 10" was used to 
find any positional coincidence bewteen the corresponding 
SN and sources in the wavdetect source list. Sources with 
such coincidences became X-ray SN candidate and went 
through another 10" box search in SIMBAD to do source 
identification. In case a candidate consists of 1 or more non- 
supernova identifications, it may be considered to be con- 
taminated by the nearby sources and would be eliminated 
from the candidate list. 

Finally, soft band (0.3-1.1 keV), median band (1.1-2.0 
keV) and hard band (2-10 keV) images were extracted from 
the event files and wavdetect with the above parameters were 
applied to each band images to check the count rates of each 
band and hardness ratios of the X-ray SN candidates. 



2.2 Justification of the searching method 

The ASDC Swift-XRT Simulatoi0 was used to simulate 
Swift observations of series of fake X-ray sources. The sim- 
ulated observations were then processed by the search- 
ing routines discussed to give blind searches of the X-ray 
sources. Percentage of sources recovered should be a key pa- 
rameter for justifying effectiveness of the program. In this 
test, an index 2 power law with hydrogen column density 
uh = 2.17 X lO^^cm"^ model was used to generate 64 X-ray 
sources in different locations in image. We simulated their 
Swift observations with exposure times 2000 seconds and 
5000 seconds. In 2000 seconds simulation, 78% of the X- 
rays sources were detected by the program while all X-ray 
sources were correctly identified and only got 3 wrong iden- 
tifications for the 5000 seconds case. 

2.3 Searching Result 

24 X-ray SNe candidates were found in the Swift public 
archive using our program in the period between Novem- 
ber 2004 and February 2011, including 14 Type II, 9 Type 
Ib/c and 1 Type la X-ray SN candidates (see Table 1). 
Many candidates have bee n previously studie d including 
SN 1 978K in NGC 1313 (|Sclilegel et al. 1 Il999l ). 1993J in 
M81 (|Chandra et al.l |2009l ) and etc. We compared our re- 
sult with the previous publications of Swift and found that 
SN 200 5ke, 2008M, 2008ax, 2008ij, 2009gj, 2010F, 2010jl, 
2010jr (|lmmler et al. 2006 , 2008 , 2009a b. 2010al lbll3ldl) were 
missing in our list. The reason was mainly believed due to 
the different approaches of the SN/source detections. The 
general X-ray SNe searching methods are targeted searches 
while a blind search method was used in this work. The 
latter could be done in an automatic way easily. In order 
to due with a 5 years database (about fifty thousand ob- 
servations) automatically, blind search approach was used 
in this project and it was supposed to be the major reason 
why there is a slightly difference between the results. In- 
deed, two of the above missing SNe, SN 2010jl and 2010jr, 
were found in the out blind search. However, they are too 
close to the host galaxies, so SN 2010jl was heavily contami- 
nated by the galaxy and SN 2010jr is very hard in spectrum 
which is likely to be something else (i.e. emission from its 
host galaxy). 



^ http:/ /www. cfa.harvard.edu/iau/lists/Supernovae. html 



^ http:/ /www. asdc.asi.it/simulator/swift/ 
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Figure 1. The XRT images of energy band 0.3-10 IceV of the "missing" Swift X-ray supernovae. The green ellipses show the sources 
detected by wavdetect in this work while the red crosses show the corresponding SN locations. 



Table 1. Complete list of supernova found in Svjift 



SN Name SN RA SN DEC OBS DATE DBS END EXP SRC RA SRC DEC SEP SIGN 

(UT) (UT) (s) (arcsec) (sigma) 



SN 1970G 


14 


03 


00, 


,83 


+54 14 32.8 


2005-01-20 


01: 


:23: 


:34 


2007-04-19 


10: 


;26: 


;56 


65481 


14 


03 


00, 


,43 


+54 14 32.1 


4, 


,1 


± 


1, 


,8 


3.50 


SN 1978K 


03 


17 


38, 


,62 


-66 33 03.4 


2006-02-03 


22: 


;56: 


:19 


2009-12-20 


02: 


;42; 


;58 


49940 


03 


17 


38, 


.66 


-66 33 02.0 


1, 


.4 


± 


0, 


.6 


70.92 


SN 1979C 


12 


22 


58, 


.58 


+15 47 52.7 


2005-11-06 


16: 


:00; 


:55 


2006-03-29 


15; 


;15; 


;58 


57636 


12 


22 


58, 


.59 


+15 47 47.9 


4, 


.8 


± 


0, 


.9 


8.60 


SN 1986L 


04 


17 


29, 


.40 


-62 47 04.0 


2005-08-10 


02: 


;45: 


:26 


2007-12-12 


16; 


;42; 


;57 


50474 


04 


17 


29, 


.19 


-62 46 58.9 


6, 


.0 


± 


1, 


.3 


14.83 


SN 1987A 


05 


35 


28, 


.01 


-69 16 11.6 


2005-01-24 


01: 


;09; 


:05 


2006-01-11 


09; 


;48; 


;56 


35728 


05 


35 


28, 


.07 


-69 16 10.8 


1, 


.1 


± 


0, 


.3 


134.90 


SN 1993J 


09 


55 


24, 


.78 


+69 01 13.7 


2005-04-21 


00: 


:47: 


:46 


2011-02-17 


14: 


;35: 


;58 


293723 


09 


55 


24, 


.64 


+69 01 13.6 


2, 


.0 


± 


0, 


.6 


31.53 


SN 1996cr 


14 


13 


10, 


.05 


-65 20 44.8 


2007-03-23 


06: 


:08: 


:28 


2009-11-16 


22: 


;54: 


;55 


18492 


14 


13 


10, 


.04 


-65 20 44.4 


0, 


.4 


± 


0, 


.5 


35.33 


SN 20031x 


16 


19 


21, 


.65 


+41 05 23.8 


2008-01-10 


20: 


;39: 


:16 


2008-01-14 


22: 


;52: 


;57 


13810 


16 


19 


21, 


.70 


+41 05 25.7 


2, 


.1 


± 


1, 


.2 


17.20 


SN 2005ip 


09 


32 


06, 


.42 


+08 26 44.4 


2007-02-14 


01: 


;57: 


:14 


2009-11-19 


05: 


;28: 


;58 


25304 


09 


32 


06, 


.44 


+08 26 44.8 


0, 


.5 


± 


0, 


.5 


32.94 


SN 2005kd 


04 


03 


16, 


.88 


+71 43 18.9 


2007-01-24 


00: 


;56: 


:52 


2008-08-21 


22: 


;44: 


;58 


18168 


04 


03 


17, 


.19 


+71 43 19.2 


4, 


.0 


± 


2, 


.0 


12.90 


SN 2005nc 


18 


32 


32, 


.58 


+26 20 22.6 


2005-05-25 


00: 


:03: 


:00 


2005-06-28 


23: 


:20: 


;57 


268557 


18 


32 


32, 


.62 


+26 20 20.6 


2, 


.1 


± 


0, 


.2 


51.77 


SN 2006aj 


03 


21 


39, 


.66 


+16 52 02.1 


2006-02-18 


05: 


:13; 


:38 


2010-12-15 


23; 


;13; 


;57 


357796 


03 


21 


39. 


.69 


+16 52 00.7 


1, 


.4 


± 


0, 


.2 


87.06 


SN 2006bp 


11 


53 


55, 


.74 


+52 21 09.4 


2006-04^10 


12: 


;51: 


:21 


2007-03-31 


22; 


;10; 


;56 


87892 


11 


53 


55. 


.84 


+52 21 13.3 


4, 


.0 


± 


1, 


.2 


3.51 


SN 2006jc 


09 


17 


20, 


.78 


+41 54 32.7 


2006-10-13 


16: 


;10: 


:14 


2011-01-04 


22; 


;02; 


;56 


202441 


09 


17 


20. 


.90 


+41 54 37.6 


4, 


.9 


± 


0, 


.4 


34.98 


SN 2006jd 


08 


02 


07, 


.43 


+00 48 31.5 


2007-11-16 


17: 


;35: 


:21 


2010-06-09 


03; 


;07; 


;56 


69347 


08 


02 


07. 


.36 


+00 48 31.7 


0, 


.7 


± 


0, 


.4 


39.68 


SN 2007od 


23 


55 


48, 


.68 


+18 24 54.8 


2007-11-05 


16: 


:13: 


:08 


2008-08-23 


19; 


;02; 


;56 


65929 


23 


55 


48. 


.69 


+18 24 52.1 


2, 


.7 


± 


1, 


.3 


4.21 


SN 2007uy 


09 


09 


35, 


.40 


+33 07 09.9 


2008-01-06 


00: 


:33: 


:07 


2009-02-21 


23: 


;01: 


;55 


334518 


09 


09 


35, 


,02 


+33 07 09.0 


1, 


.1 


± 


0, 


.7 


5.82 


SN 2008bo 


18 


19 


54, 


.41 


+74 34 21.0 


2008-04-04 


14: 


:54: 


:02 


2009-02-14 


23: 


;05: 


;57 


383873 


18 


19 


54, 


,60 


+74 34 20.7 


2, 


.0 


± 


0, 


.6 


34.77 


SN 2008D 


09 


09 


30, 


.62 


+33 08 20.1 


2008-01-06 


00: 


:33: 


:07 


2009-02-21 


23: 


;01: 


;55 


334518 


09 


09 


30, 


,62 


+33 08 20.2 


0, 


.0 


± 


0, 


.3 


46.80 


SN 2008hw 


22 


39 


50, 


.39 


-40 08 49.1 


2008-10-07 


05: 


:26: 


:52 


2008-10-27 


14: 


;00: 


;56 


188541 


22 


39 


50, 


,59 


-40 08 47.0 


3, 


.1 


± 


0, 


.3 


85.51 


SN 2009dd 


12 


05 


34, 


.10 


+50 32 18.6 


2009-04-15 


14: 


:50: 


:05 


2009-09-21 


23: 


;59: 


;57 


104014 


12 


05 


34, 


,58 


+50 32 20.7 


7, 


.3 


± 


1, 


.0 


8.85 


SN 2009mk 


00 


06 


21, 


.37 


-41 28 59.8 


2009-12-17 


03: 


;30; 


:58 


2010-01-30 


23; 


;59; 


;57 


38301 


00 


06 


20. 


.72 


-41 29 03.8 


8, 


.8 


± 


1, 


.4 


1.90 


SN 2009nz 


02 


26 


19, 


.90 


-18 57 09.0 


2009-11-28 


01: 


;02: 


:16 


2010-06-10 


23; 


;44; 


;57 


470379 


02 


26 


20. 


.03 


-18 57 07.2 


2, 


.6 


± 


0, 


.1 


154.47 


SN 2010bh 


07 


10 


31, 


.80 


-56 15 20.2 


2010-03-16 


21: 


;53: 


:21 


2010-03-23 


20; 


;55; 


;33 


52488 


07 


10 


30. 


.95 


-56 15 17.5 


12.5 


i ± 1.2 


13.18 
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2.4 Data Reduction 

We used HEAsoft of version 6.6.2 to do the data analy- 
sis. Tasks xrtgrhlc and xrtgrblcspec were used to extract 
lightcurves and spectra while the Swift-'KKT data products 
generator (j Eva ns ct al. 2009) was use^fl to do some cross- 
checks sometimes. We binned all spectral products with 20 
counts per bin at first. Those binned spectra with 4 or less 
data points were re-binned with fO counts per bin. Re- 
binned spectra with 7 or less data points were re-binned 
with 5 counts per bin. Spectra with 9 or less data points 
after 5 counts binning were regarded as low signal to noise 
which were not used to do model fitting. The background 
regions were chosen automatically by the tasks. Source-free 
background regions would be extracted by human some- 
times to avoid contamination from some nearby bright X-ray 
sources or inappropriate background selections at the edges 
of images. Finally, XSPEC was used to do the spectral fit- 
ting to examine physical parameters of the sources. Basi- 
cally, combinations of simple photon p ower law model and 
Raym ond-Smith thermal plasma model (|Ravmond fc SmithI 
1 19771 ) were primary models of the analysis. Column densi- 
ties from Leiden/ Argentine/Bonn (LAB) Survey of Galactic 
HI were used as the Galactic uh, and uh excess components 
would be added if necessary. The best fit results of chi-square 
fittings of 16 X-ray SN candidates are shown in the Table [2] 



2.5 Properties of X-ray supernovae spectra 

In the early phase of the explosion, the strong photosphere 
radiation is Compton scattered which forms a high energy 
tail to the photosphere emission. Each photon increases its 
energy by a factor Au/iy = AkTc/iricC^ > 1 where nie and 
Te are electron mass and electron gas temperature, which 
create a power law spectrum with a photon index ((Franssonl 
ll982l : IChevalier fc Franssoiill200ll ') 



Q = {- 



ln[y (0.9228 



Inre)]}^ 



(1) 



where Te is the optical depth to electron scattering behind 
the circurnstella r shock. The typical value is 1 < a < 3 
iFranssonl Il982l '). In order to compare the SNe data to 
the model prediction, hardness ratio diagram of high sig- 
nal (count > 50) X-ray supernovae candidates was plotted 
(Table [2|. In the diagram, hardness ratios HRi — j^~g and 
HR2 = were plotted where S = soft band (S: 0.3-1.0 
keV), M = median band (M: 1.0-2.0 keV), and H = hard 
band (H: 2.0-10 keV). With the reference trends from the 
Portable Interactive Multi-Mission Simulator (PIMMS) pre- 
diction of power law model with various photon indexes and 
column densities, the diagram shows clearly that the ratio 
between the hardness HRI and HR2 are consistent with a 
power law model of index ~ 2 which is consistent the predic- 
tion. Spectral fitting resuls also agreed with the theoretical 
model. As shown in the Table [S] most of the supernovae (ex- 
cept SN 1993J) show non-thermal components (power law) 
in their spectra with a mean photon index ~ 2 while some 
spectra also show thermal components with plasma tem- 
peratures about ~ 0.8 keV. If ages of SNe were taken into 



http:/ /swift. ac.uk/user_objects 



Swift Data 
PL with index 1 .5 
PL witti index 2.0 
PL with index 2.5 



HR1 [(IVI-S)/(IVI+S)] 

Figure 2. Hardness diagram of the X-ray SNe. The points are 
the PIMMS of version 4.1a predicted values using Power Law 
(PL) with different values of column density (from left: 0.1, 1, 2, 
3, 4xl02lcm-2). The soft band (S: 0.3-1.0 keV), median band 
(M: 1.0-2.0 keV), and hard band (H: 2.0-10 keV) hardness ratios 
HRi and HR2 are plotted. 



the account, it is clear that older supernovae spectra were 
fit well with thermal plus non-thermal model while younger 
supernovae could be satisfied with single power law. One 
explanation is that non-thermal blast shock emissions dom- 
inated the whole X-ray emissions of supernovae which cre- 
ated single power law continuum in the early phase. After 
years, the shock velocities decreased as the swept-up mass 
behind the circumstellar shock increased which weakened 
the non-thermal part, hence the thermal emissions (and the 
line emissions from elements) started to be significant. 



2.6 X-ray Supernovae luminosity 

We estimated the X-ray luminosity of the X-ray super- 
nova candidates by the spectral model fitting. For those low 
signal-to-noise candidates data, the net count rates with 
power law of photon index 2 were used to estimate the 
flux using PIMMS. Fig. [3] shows the luminosity function 
of all candidates. Most of candidates have a luminosity of 
~ lO^^erg/s in energy band 0.3-10 keV. The most luminous 
SNe were the GRB-SN associations which have a luminos- 
ity of ~ 10''*erg/s while the dimmest one was ~ 10='''erg/s, 
which belongs to the SNR 1987A in Large Magellanic Cloud. 
This wide X-ray luminosity range partly reflects the different 
stages of evolution of the observed supernovae. 

In addition to know more the properties of the SNe 
progenitors, we used the X-ray luminositie s to es timate the 
CSM density profile. [Chevalier fc FranssonI l|l994h suggested 
free-free luminosity from reverse shock can be estimated 
from 



Lrcv 



ttR^ nc{R) ni{R)A{T)tV 



■ ergs 



(2) 



l + [2nc(R)A{T)/{mp V^)]t 

where R and T are radius and temperature of the reverse 
shock, ne and rn are the number densities of electrons and 
ions, t is age of the supernova, V is the shock velocity, rrip 
is the proton mass and A is the cooling factor, which can be 
expressed as A = 6.2 x 10"^^ T'" '^ ergcm^ s"^ for 10^ K < 
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Table 2. Spectral properties of X-ray SNe candidates 



SN Name 


Distance 


njj (Galactic) 


njj (Intrinsic) 


Photon Index 


Plasma Temperature (kT) 


Reduced-x^ 




(Mpc) 


(lO^O cm-2) 


(10^0 cm-2) 




(keV) 




SN 1978K 


4.5 


4.07 


24+^ 


2 7+?,1 
^' ' -0.4 


"■^—0.1 


1.08 


SN 1986L 


16 


2.17 




— 0.2 


0.7+11 

— V.I 


1.25 


SN 1987A 


0.05 


26.5 




3 s+V 


7+5? 

" ' -0.1 


1.42 


SN 1993J 


3.6 


5.57 






0.8+0-; 

— U. 1 


1.04 


SN 1996cr 


17 


55.6 




1 1+"-2 

2 




1.27 


SN 2005ip 


30 


3.56 




5+°'^ 




1.16 


SN 2005nc/GRB050525A'' 


1780 


9.07 




1 7+0.1 
^- '-0.1 




1.04 


SN 2006aj/GRB060218'' 


136 


9.37 


46+"^ 
™-10 


n + 0.4 
■^■°-0.3 


7+°-i 

'^•'-0.2 


1.58 


SN 2006jc 


23.6 


1.51 


19+13 


1 7+0.2 
^■'-0.2 




0.81 


SN 2006jd 


79 


4.52 




1 9+0.2 




1.33 


SN 2008bo 


21 


5.37 




1 9+°-2 




1.16 


SN 2008D" 


27 
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a. Spectrum after outburst was used. 

b. Time average spectra were used. 



T ^ 4 X 10^ K. With the assumption of the CSM is hydrogen 
and R2 = Vc]t = {n — 2) V2 t, we get 



2/iH2me2(n-2){l + 



4AprcvMc (n-2)^t 



ergs 



(3) 



where pcs and p, eu refers to the densities behind the circum- 
stellar shock and the reverse shock in unit of g cm~'^ . We 
assumed the progenitors are Wolf-Rayet stars with a mass 
loss rate M = 2 x 10"'*MQyr"i and w„ind = lOOOkms"^ 
Some typical values n — 15 and T = 10^ K (which is about 
0.7 keV which is consistent to our spectral fitting result) 
and Vcj = 35000 km s~i were also applied in the calcula- 
tion. With the relation prcv/pcs ~ 0.602 (n — 3)^ and the 
X-ray luminosities of various Type Ib/c SNe, we get a typ- 
ical picture of Type Ib/c CSM density profile. After com- 
paring it with a steady flo w model, pcs = M /nw^indr'^ 
(jChevalier fc FranssonI [l994l ) . the difference between them 
increases as time/radius raises. If a dropping temperature 
(i.e. T oc t~^'^) was used , the result then matched perfectly 
with the prediction of the steady flow model. Although the 
calculation looks perfect, it is important to notice that the 
decay rate T oc t~^'^ is not signiflcant to the real situation 
since there are many assumptions and fixed parameters (i.e. 
the mass loss and the shock velocity) in the model. The de- 
cay rate would be different if other parameters were freed. 



3 STUDY OF INDIVIDUAL X-RAY 
SUPERNOVA 

3.1 SN 1986L 

SN 1986L was discovered on October 7.6, 1986, with an 
apparent magnitude of 13.7 in the location of 50" due 
west of the NGC1559's nuc leus (R A. = 4hl7m0, Peel. = 
-62deg55', equinox 1950.0) (|Evanslll986l : lMcNaughtlll986f ). 
It is located at NGC1559 where is 16 Mpc from the earth 




38 40 42 

LOG[Luminosity„,,„,J 

Figure 3. The observed X-ray 0.3-10 keV luminosities of su- 
pernovae calculated by PIMMS, from SN 1987A (blue shadow) 
at logio(Lo.3-iokev) = 36.3 to SN-GRB candidate SN 2008hw 
at log^Q(Lo.3— lokev) = 44.2, with an average peak luminosity 
l°gio(Lo.3— lOkov) = 40.3, reflecting partly the different stages of 
evolution of the observed supernovae. 



(|Schmidt fc Kirshner et~allll994l ). SN 1986L was classified 
as a core-collapse Type II supernova. 

Swift did several observations on the field of SN 1986L 
from 2005-08-10 to 2007-12-12 and detected an excess of X- 
ray counts (0.3-10 keV) located on 7.1" from the optical 
SN position at a 15.4 a level in the 50.5 ks combined image. 
The XRT net count rate is (2.5 ± 0.3) x 10"^ count/s. By 
assuming the column density uh = 2.17 x lO^^cm"^ and its 
distance from the earth was Mpc, the luminosity was ap- 
proximately (1.8lo;4) X 10*°ergs"i in 0.3-10 keV band. 

The ROSAT High Resolution Imager (HRI) and XMM- 
Newton observed at SN 1986L on 1990-06-01 and 2007-04-26 
with exposure time 17.5ks and 12.5ks, respectively (Fig. |4]). 
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Figure 4. Aligned soft (0.3-2.4 keV) X-ray images of SN 1986L taken by ROSAT (Left), Swift (Middle) and XMM-Newton (Right) with 
the position of the supernova indicated by a 10 arcsec radius green circle. The sUghtly different positions of the X-ray source observed 
are consistent with the astrometric accuracies of the different satellites. 



In the figure, A position coincidental X-ray source was ob- 
served in all ROSAT, Swift data and XMM-Newton data. 
The sources seem to have different locations in each mis- 
sion due to the different position accuracies (i.e. ~ 6" for 
ROSAT- ~ 5" for Swift; ~ 2" for XMM-Newton). The low 
exposure time of ROSAT observation brought a low signal- 
to-noise ratio image which is not much informative. Using 
the data of XMM-Newton, the luminosity in energy band 
0.3-10 keV was estimated which is (S.ll?;^) x 10^® erg s"^ 
in 2007. 



3.2 SN 1987A 

SN 1987A is one of the most famous and well-known su- 
pernovae in astronomical history because it is the only su- 
pernova that was detected in neutrino before optical light 
reached the earth till now. SN 1987A is located in the Large 
Magellanic Cloud, which is 50 kpc from us. This short dis- 
tance allows us to resolve the detailed structure of the en- 
tire debris in optical. Apart from the optical features, the 
SN 1987A is distinctive because it is the only Type IIP su- 
pernova that was detected in X-rays. It is also the X-ray 
dimmest supernova (~ 10'^^ erg s~^ in 0.3-lOkeV) ever de- 
tected, which is the result of the sh ort distance between 
LMC and the earth. (|Triimpeij[200i ). The soft band (0.5- 
2.0 keV) light curve of the SN 1987A in the past ~10 years 
(Fig. [SJ shows that an exponential increasing trend almost 
from ~5500 days after the outburst. This rapid increase was 
due to the interaction between the blast shock wave and the 
medium in the inner ring. 

Swift observed SN 1987A six times between April 
2005 and January 2006. Using a double constant 
temperature plane-parallel shock plasma model (i.e. 
wabs*vpshock*vpshock) with hydrogen column density 
and elementa l abundances from the best fit values of 
IZhekov et al.l (|2006l ). the best fit shock temperatures are 

keV respectively which are 



,+0.10 
-0.10 



7+0.60 
-0.74 



consist with results 0.411q25 



keV and 3.40to.66 keV of 



jZhekov et al.ll2006l ). The estimated luminosity of 0.5-2.0 
keV is L - 1 
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erg s ^ which is also consistent 



Figure 5. X-ray light curve of SN 1987A in the 0.5-2.0 keV band 
for the past ~ 10 y ears. The XMM-Newton flux val ues are from 
iHab erl et al.l l|2006h (for the first three data) and iHeng et al.l 
(i2008l) (for the l ast data). The Chandra flux values are from 
iPark et al. I | |2007|) . The Swift fluxes are from this work, which 
were estimated by the XSPEC model wabs*(vpshock-(-vpshock) 
with temperatures of 0.3 and 2.3 keV, and abundances from 
IZhekov et all bOOdt) . 



to the trend of the multi-mission lightcurve of SN 1987A 
which is shown as Fig. (5) 



3.3 SN 1993J 

SN 1993J was disco vered on March 28, 19 93 in the nearby 
spiral galaxy M81 (|Ripero &: Garci i|l993'). It is the most 
intensively monitored X-ray SN after SN 1987A. The first X- 
ray observation was done by ASCA 6 days after its discovery. 
Afterward, the evolution of the SN was monitored by ASCA, 
ROSAT (6-1800 days); XMM-Newton, Chandra (2500-5500 
days); Swift (after 4800 days) in the past 16 years. With 
the substantial number of complete observations, the entire 
SN evolution analysis could be performed through spectral 
fitting. 
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a rapid enhancement of X-ray flux. As shown in the Fig. [6l 
the X-ray luminosity of the SN 1993J was decreasing over its 
hfe. IZimmermann et al.l (|l993l ) suggested the luminosity in 
range of 0.3-2.4 keV decay with t~'^ over the first 200 days. 
With the assumptions of the mass loss rate of the progenitor 
wind can be expressed as 



(4) 



and the luminosity is roughly proportional to square of the 
density, the density profile index w as estimated as s = 1.65 
IZimmermann fc Aschenbachll2003l ). 

The Swift data taken between 2009 and 2011 (Fig. ^ 
shows the decay rate of the X-ray flux was slowing down. 
The light curve is almost a constant ??? or ??? with one 
sigma fluctuation ??? or ??? in that period. 



Figure 6. 0.3-2.4 keV X-ray lightcurve compiled from ASCA: tri- 
angles; ROSAT; squares and pentangles; Chandra: circles; XMM- 
Newton: crosses; Sw ift: x crosses are original data reduced by 
I Chandra et al Swift in 2009: filled circles are original in 

this work. 
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Figure 7. Swift X-ray lightcurve of energy band 0.3-10 keV of SN 
1993J in 2009. The curve is flat and the count rate is low compare 
with the former years. There was a short bump in the X-ray light 
curve after day 6025 which shows the uneven distribution of the 
CSM. 

The early detection suggested that the existence of 
dense CSM region which was close to progenitor. OSSE 
on the GRO satellite determined the init ial temperature 
of SN 1993J on the order of 80 keV (|Triimpeil l2008l : 
IZimmermann fc Aschenbachll2003l ) in the flrst month. With 
the extreme high initial temperature and the limitations of 
the instrument, the ASCA meeisurements a lower lim it of 
10 keV was reported (jZimmermann fc Aschenbachll2003f ) . In 
April 2001, SN 1993 J was observed with both the PN camera 
of the XMM-Newton observatory for exposing about 70 ks. 
The XMM-Newton spectrum was fltted by Ao\ih\e-vmekal, 
a 2-component thermal model. The best fit indicated that 
the shock wa ve temperatures were 0.34ln'n^ and 6. 541* keV 
respectively l|Zimmermann fc Aschenbachl200^ ) after about 
2900 days. 

Unlike the unusual X-ray evolution of SN 1987A with 



3.4 SN 2007od 

SN 2007od was detected by S. Maticic on four unfiltered 
CCD images taken around Nov. 2.85 UT with the 60-cm 
f/3.3 Cichocki reflector in the course of the Comet and As- 
teroid Search Program (PIKA) at Crni Vrh Observatory. 
The optical posi tion is at R.A. = 23h5 5m48s.68, Deck = 
+ 18deg24'54.8" (|Mikuz fc Maticicll2007l ). 

SN 2007od was observed by Swift initially 3 days af- 
ter the explosion. No significant X-ray signal was detect at 
the position in the 4.0 ks XRT exposure by that observation 
(|lmmler et al.ll2007b, ') . After combining with the later obser- 
vations, there was a X-ray source detected at the SN position 
with a 5.7(T significance in the merged 58.2 ks XRT data. The 
detected count rate was (4.8 ± 1.2) X lO"'' count s"^ which 
is corresponding to a flux rate of 1.8 x 10~^* erg cm~^s~^ 
with a power law model of photon index 2. The reason of 
the deferred X-ray emission is most likely absorption by the 
cool, d ense shell formed by the reverse shock at the early 
epoch dChevalier fc Franssoniil994 : iNvmark et al.ii2006 ). 

Lx OC J-(15-6s-|-sn-2n)/(ns) j-g-j 

where r~" is the radial dependence density profile of the SN 
ejecta and s is the profile index of (4). By taking a reason- 
able assumption of th e progenitor is compact with n = 10 
l|Fransson et all 1 19961 ). it gives oc t-(4«-5)/(io<i) 
pare the equation with the observed decline of t~'^'^ , the 
value s ~ 2 which is consistent to the initial interaction of 
an expanding Type II sup ernova with a CS M created by 
mass loss of the progenitor (|Chevalieilll982bl ). 



4 DISCUSSIONS AND CONCLUSIONS 

We have investigated over 6 years archival Swift observations 
taken from November 2004 to February 2011 and discovered 
three new candidates (i.e. SN 1986L, 20031x and 2007od) out 
of 24 X-ray detectable Swift SNe. It is considerable because 
there were only ~ 40 detections of X-ray SNe in the past, 
we recovered half of them and increased the total detections 
by ~ 10%. Furthermore, The automatic selection program 
developed in this project worked well to select useful su- 
pernova data from over thirty thousand observations in the 
archive. The success of this work ensures the feasibility of 
other analogous surveys in Swift. It is also promising that 
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Figure 8. Swift X-ray lightcurve of SN 2007od in range of 0.3 to 
10 keV. 



more X-ray SNe detections would be found if this complete 
survey lasts for a few more years. 

There is no doubt that every supernovae has its unique 
features. Even the supernovae of the same type, the varia- 
tions between environments, stellar abundances, stellar wind 
and shock wave properties are doubtless to influence the X- 
ray emissions. It is difficult to find a simply model to satisfy 
all features of all X-ray SNe spectra perfectly. However, in 
a macroscopic view, it is still fair to say that most of the X- 
ray SNe spectra are similar to each other. Power law model 
matches most of the cases may be a tangible testimony of 
that. Based on this thought, we have done several statistics 
on their luminosities, hardness ratios and the CSM density 
profile. The results show that most of the SN luminosities 
lie on the 10^^ to lO^^ergs"^. The studies of hardness ra- 
tios and CSM density profile partly ensure the the oretical 
models (|Franssonlll9^ : IChevalier fc Franssonlll994 l. 

One of the important contributions of this work is the 
discovery of the Type la X-ray SN candidate, SN 20031x. 
The X-ray luminosity was ~ 4.8 x 10*^ erg s"^ in 0.3-2.0 
keV which was inferred by power law of index ~ 2 and 
distance 155 Mpc. The luminosity is high compare to the 
other la SNe whic h have typical luminosity of a upper limit 
10 38-39 ijimmler et al.l l2006: Hu ghes fc John et all 

l2007h . However, observations taken by ROSAT in 1990 pro- 
vides a strong evidence that there was no X-ray detection at 
the SN position before the SN explosion which is consistent 
to a Type la SN emission. 

We did some study on the SN 1987A, 1993Jx to monitor 
the their evolutions. Basic analysis of the three new SNe 
were done to see their X-ray properties. In fact, there are 
still several X-ray SNe potential candidates (for example, 
SN 1989B and SN 2007ax) which did not show in this paper 
as their X-ray signals are weak and/or the off-set from SN 
position are huge. More Observations or better images (by 
Chandra or XMM-Newton) are expected to confirm these 
potential candidates. 
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